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We summarize recent progress in soft QCD modeling based on the set of Dyson— 
Schwinger equations truncated to ladder-rainbow level. We pay particular atten- 
tion to electromagnetic elastic and transition form factors of the pion. This covari- 
ant approach accommodates quark confinement and implements the QCD one-loop 
renormalization group behavior. The dressed quark propagators are compared to 
the most recent lattice-QCD data. 



1 Introduction 

High energy electroweak probes utilize perturbative QCD to quantify the 
structure of hadrons in terms of parton structure functions. Their detailed 
behavior samples various aspects the quark-gluon dynamics of the hadron 
bound state. This information is intrinsically non-perturbative and the un- 
raveling of it from sets of structure functions will require a connection be 
made with non-perturbative QCD calculations and models. Although sonoe 
lattice QCD studies have begun to produce moments of structure functions 
the opportunities are very limited at present. The most extensive hadronic 
models were designed to study the mass spectrum and decays and usually con- 
tain elements that limit their use for high energy lepton scattering. Examples 
include: non-relativistic kinematics, a lack of manifest Lorentz covariance, no 
quark sea, no dynamical gluons, no renormalization group behavior for change 
of scale, and no confinement of quarks. Studies of deep inelastic scattering 
within the Nambu-Jona-Lasinio model □ have helped clarify some of the issues 
confronting work within a quark field theory format. 

Here we summarize recent progress in soft QCD rnDdeling based on the 
set of Dyson-Schwinger equations [DSEs] of the theory Wa; we pay particular 
attention to electromagnetic couplings to the pion. This covariant approach 
accommodates quark confinement and implements the QCD one-loop renor- 
malization group behavior. Deep inelastic scattering phenomena have not 
been treated with the present model; a simplified related approach that re- 
tains all the essential features has recently jsroduced excellent results for the 
pion valence quark distribution amplitude u. 
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Meson electromagnetic form factors in impulse approximation are de- 
scribed by two diagrams, in one the photon couples to the quark, in the 
other to the antiquark. A form factor can be associated with each diagram. 
For example, thcjihoton coupling to the antiquark of a flavor ab pseudoscalar 
meson produces □ 

2P^F,-,-M') = N,l''^TT[S^iq)TfMq+:P^)S'{q+) 

XiTl{q+,q^)S\q^)r^liq^,q;-P+)] , (1) 

where q = k+\P, q± = k-\P±\Q, P± = P ± \Q. We work in the 
isospin symmetry limit, so for the pion we have F^{Q'^) — F^^uuiQ"^)- The 
charged and neutral kaon form factors are given by Fx+ = f^lisu + ^Fuss 
and Fj(o = —^Fdgd + j^dss respectively. To implement such a study, one 
needs a consistent QCD modeling of the photon-quark dressed vertex T^, the 
dressed quark propagator S{q), and the meson Bethe-Salpeter (BS) amplitude 
Tps. The first two are related by the vector Ward-Takahashi identity [WTI]; 
the second two are related by the axial-vector WTI. 

The DSE for the renormalized quark propagator in Euclidean space is 

S{p)-' =tZ2^ + Zimif,) + Zi j^-^g^D^,{k)^^^S{q)Tl{q,p), (2) 

where D^,^{k) is the dressed-gluon propagator, rj5(<7;p) is the dressed- 
quark-gluon vertex, and k = p — q. The solution has the form 
S{p)~^ — i-^A{p^) + B{p^) and is renormalized at spacelike according to 
= 1, and B{^^) = m(/Lt), where to(/x) is the current quark mass. The no- 
tation represents a translationally invariant regularization at scale A. The 
same regularization occurs consistently at all stages as indicated in Eq. (|l|). 
One takes A — > oo as the final step. 

The qq^ vertex r^(p+,p_) corresponding to total momentum 
Q = p+ — p^ satisfies the Bethe-Salpeter equation [BSE] 

where we use the notation p+ — P + r]Q and p- = p — (1 — i])Q foi" the out- 
going and incoming quark momenta respectively at the vertex. Here 77, which 
specifies how the total momentum is shared between quark and antiquark 
(thus defining a relative momentum p) , is arbitrary and physical observables 
should not depend on it. The kernel K is the renormalized, amputated qq 
scattering kernel that is irreducible with respect to a pair of qq lines. Other 



r^b+iP-) = ^2 7m + / K{p,q;Q) S{q+)T^{q+,q-) S{q-) , (3) 
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electroweak processes Q require, e.g., the dressed qqW vertex which is defined 
in an analogous way. 

When the homogeneous version of Eq. (|^) has vector solutions at (dis- 
crete) timelike total momenta — —m?, these are the vector meson bound 
states. The vector vertex r^(p+,p_) will have simple poles at those lo- 
cations. This correspondance holds for all other transformation characters 
(pseudoscalar, axial vector, etc), labelled by quantum numbers J^'~^ . The 
homogeneous solutions T m{p+tP-) are the bound state BS amplitudes and 
they are normalized in the canonical way. Explicit representations of the 
T m{p+tP~) require an expansion in a complete set of covariants constructed 
from gamma matrices and momenta. For example, the general representation 
for pseudoscalar bound states is eI 

Tps{k+,k^:P)^-t^[iE+ fF+ ftG + a^,k^P,H] , (4) 

where the invariant amplitudes E, F, G and H are Lorentz scalar functions 
f{k^; k ■ P; if). Note that these functions depend on the momentum partition- 
ing parameter rj because the total and relative momenta have been employed 
in the covariants. However, physical observables are independent of this pa- 
rameter; this is verified numerically within, the present approach as long as 
the set of employed covariants is complete u. 

2 Ladder-Rainbow Model 

To solve the BSE, we use a ladder truncation, with an effective quark- 
antiquark interaction that reduces to the perturbative running coupling at 
large momenta EfLI. In conjunction with the rainbow truncation for the quark 
DSE, the ladder truncation of the BSE preserves both the vector WTI for the 
gg7 vertex and the axial-vector WTI. The latter ensures the existence of al- 
most massiess pseudoscalar mesons which are the Goldstone bosons connected 
with dynamical chiral symmetry breaking u. In combination with the impulse 
approximation, this ladder-rainbow truncation ensures electromagnetic cur- 
rent conservation 

The ladder truncation of the BSE, Eq. (|), is K{p,q;P)^ 

-a(fc2)-Dj;r(fc) T^A'^T^ ^here k = p - q, and D^-'^'' is the free gluon 
propagator in Landau gauge. The corresponding rainbow truncation of 
the quark DSE, Eq. (||), is given by T^^{q,p) ^ /2 together with 
9^Df,^{k) g{k'^)D^ff^''{k). This truncation was found to be particularly suit- 
able for the flavor octet pseudoscalar and vector mesons since the next-order 
contributions in a quark-gluon skeleton graph expansion, have a significant 
amount of cancellation between repulsive and attractive corrections E2I. 
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Table 1. Results for light mesons in GeV compared to experiment 





expt. (estimates) 


calc. ( fitted') 


''*^=lGcV 


P' 1 /~\ TV T "Y y 

5-10 MeV 


5.5 MeV 




100 - 300 MeV 


125 MeV 


- imr. 


(0.236 GeV)^ 


( 0.241 GeV)^ 




0.1385 GeV 


0.138 


U 


0.0924 GeV 


0.093 


rriK 


0.496 GeV 


0.497 


fK 


0.113 GeV 


0.109 


Wlp 


0.770 GeV 


0.742 


fp 


0.216 GeV 


0.207 


rriK* 


0.892 GeV 


0.936 


fK* 


0.225 GeV 


0.241 


f<P 


1.020 GeV 


1.072 


0.236 GeV 


0.259 




0.44 fm^ 


0.45 


^K+ 


0.34 fm^ 


0.38 




-0.054 fm^ 


-0.086 


' 7777 


0.42 ivc? 


0.39 


^7777 


0.50 


0.50 




0.57 


0.54 




6.02 


4.85 


g<pKK 


4.64 


4.63 


9K* + tt+K» 


4.53 


4.6 



The model is completely specified onne a form is chosen for the "effective 
coupling" G{k^ )- We employ the Ansatz l3 



fc2 



47r2 D e 



47r2 7„^(fc2) 



In 



+ (l + fcVA2jcD) 



(5) 



with 7„ = 12/(33 - 2Nf ) and T{s) = (1 - exp ^)/s. The ultraviolet be- 

havior is chosen to be that of the QCD running coupling a(A;2); the ladder- 
rainbow truncation then generates the correct perturbative QCD structure 
of the DSE-BSE system of equations. The first term implements the strong 
infrared enhancement in the region < fc^ < 1 GeV^ phenomenologically re- 
quired to produce a realistic value for the chiral condensate. In recent years, 
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there has been some progress in understanding the interplay between ghosts 
and gluons in Landau-gauge QCD, indicating that ghosts could play an im- 
portant role for this enhancement □. However, this is not yet ready for 
incorporation into phenomenological studies of hadron properties. We use 
nit = 0.5 GcV, r = e^ - 1, Nf = 4, Aqcd = 0.234 GeV, and a renarmaliza- 
tion scale fi — 19 GeV which is well into the perturbative domain EnJ. The 
remaining parameters, uj = 0.4 GeV and D — 0.93 GeV^ along with the quark 
masses, are fitted to give a good description of {qq), m^/x and The value 
of Jk, the ground state vector masses and decay constants, and the ^fiu^ 
quantities in Table |^ are then produced without any free parameters Qu^M. 
Note that the complete set of covariants of the BS amplitudes, see Eq. 
are needed to satisfy the axial- vector WTI U. The pseudovector amplitudes F 
and G contribute about 30% of the masses and decay constants of the pseu- 
doscalaxs; they dominate the asymptotic behavior of the pion charge form 
factor Ej. 




Figure 1. Lefh Panel: DSE solution u for quark propagator amplitudes compared to recent 
lattice data tj. RigM-iPanel: The chiral limit DSE mass function compared to the lattice 
chiral extrapolation 



In Fig. m we compare the DSE model □ propagator amplitudes defined Jay 
S{p) = Z{p^)[i^ + M(p2)]-i with the mflst recent results from lattice QCD 13. 
In contrast to a previous comparison £3, the data here has more lattice arti- 
facts removed and the DSE calculations use the indicated current mass values 
to match the lattice mass function at 6 GeV. There is agreement in the qual- 
itative infrared structure of the mass function particularly in the way the 
infrared enhancement sets in. Since the lattice simulation produces the reg- 
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ulated but un-renormalized propagator, the scale of the field renomialization 
function Z is arbitrary and only the shape is a meaningful comparison. For 
this reason, we have scaled the lattice data for Z so that Z(5 GeV) — 1. The 
ladder-rainbow DSE model typically produces a Z that saturates much slower 
than does the lattice Z; tkis may signal a deficiency of the bare gluon-quark 
vertex. The lattice work IIj also produced a linear extrapolation to the chiral 
limit mass function Mo(p) and we compare this to the DSE result with the 
insert emphasizing high momentum. The known one-loop renormalization 
group UV behavior of chiral QCD, which is preserved by this DSE model, is 



27r27„ 



In 



TP — 



l-7r. 



(6) 



with {qq)'^ being the renormalization-point-independent chiral condensate la. 
Although the lattice data does not distinguish clearly between this fall-off and 
a flat behavior, the scale of the data is consistent with a condensate that is less 
than 10% higher than that of the DSE model. That is, we make a tentative 
assignment of {qq)^^i gcV = —(248 GeV)^ for this lattice data. 

3 Meson Elastic and Transition Form Factors 

Our results i for i^^((5^) are shown in Fig. H the charge radii are given in Ta- 
ble|l|. Up to about = 3 GeV^, our results for F^{Q'^) and Fk{Q'^) can be fit- 
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(GeV/c)^ 

Figure 2. Our results for the pion charge form factor ^ compared to new data from JLab 
including a re-analysis of older data. Aleoi shown are various earlier model results. Figure 
taken from the PhD thesis of J. Volmcr 
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ted quite well by a monopoleH. Asymptotically, the behavior is Q^F{Q^) — > c 
up to logarithmic corrections. However, numerical limitations prevent us from 
accurately determining these constants. Around — 3GeV^, our result for 
Q^F^ is well above the pQCD result B Wnf^asiQ^) ~ 0.2 GeV^ and clearly 
not yet asymptotic. The calculated time-like electromagnetic form factor dis- 
plays vector meson bound state poles from which_pnc can extract the strong 
coupling constants ^ptttt and g^KK respectively t3. The results are given in 
Table I ancLare close to the experimental data and to the results from direct 
evaluation Ej. 

The impulse approximation for the 7*7r7 vertex with 7* momentum Q is 

= T j J^i'^^[S{q)lT,{q,q') S{q')lT^{q' ,q") S{q")T.{q'\q;P)] . 

where the momenta follow from momentum conservation. In the chiral limit, 
the value at = 0, corresponding to the decay 7r° — > 77, is given by the 
axial anomaly and its value = i is a direct consequence of only gauge 
invariance and chiral symmetry; this value is reproduced by our calculations 113 
and correspondsjsjell with the experimental width of 7.7 eV. In Fig. || (left) we 
show our results 113 with realistic quark masses, normalized to the experimental 
f/7r77- For such moderate Q^, both the data and our DSE model are close to the 




Figure 3. Left Panel: The 7* 777 form factor B, with data from CLEO and CELLO E|. 
Right Panel: The form factor for 7* 7r7* for equal virtuality with the insert showing the 
asymptotic behavior of Q^F{Q^). 
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monopole shape fitted to the Brodsky-Lepage asymptotic Hmit 
from pQCD in the factorized approximation. As pointed out recently e3, it 
is advantageous to consider the process wher^L-bcith photons are virtual; the 
light-cone operator product expansion yields cd'Ea 

Here the asymmetry is a; = {Ql — Q2)/{Qi + Qi); ^^'^ the leading coefficient 

jH^-y^ ^^ i_,2(2,„i)2 ^ (9) 

samples the pion bound state through the distribution amplitude 4>Tr (x) . The 
normalization of the latter produces the model-independent result J(0) = 4/3 
for the case of equal virtuality. The model-dependent J(l\_£pr the asymmetric 
case has received the estimates: Jl from use of 0^^'""(a;) t3, 1.8 from a non- 
perturbative chiral quark model Ej, 1.6 from^£± to CLEO data c3, and 4/3 
from analysis of the impulse approximation c3 c3. A numerical confirmation 
of the asymptotic behavior within the present DSE model requires a level of 
detail in the representation of the analytic structure of dressed propagators 
and BS amplitudes that is not currently available. However, for equal virtu- 
ality there is no restriction, and our asymptotic result shown in Fig. ^ (right) 
produces J(0) = 1.2 which is within 10% of the 4/3 result of Eq. from the 
light-cone operator product expansion. 
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